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ABSTRACT 

The observed abrupt torque reversals in X-ray pulsars, 4U 1626-67, GX 1+4, and 
OAO 1657-415, can be explained by transition in accretion flow rotation from Keplerian 
to sub-Keplerian, which takes place at a critical accretion rate, ~ 10^^ — 10^'^ g/s. 
When a pulsar system spins up near equilibrium spin before the transition, the system 
goes into spin-down after transition to sub-Keplerian. If a system is well into the 
spin-up regime, the transition can cause a sharp decrease in spin-up rate but not a 
sudden spin-down. These observable types of abrupt torque change are distinguished 
from the smooth torque variation caused by change of M in the Keplerian flow. 
The observed abrupt torque reversal is expected when the pulsar magnetic field 
i?* ~ 5 X 10^^6p ^^"^ where the magnetic pitch parameter hp ^ a. few, Lx^i,& is 
the X-ray luminosity in IQ^^erg / s, and -P*,io is the pulsar spin period in 10s. Observed 
quasi-periodic oscillation (QPO) periods tightly constrain the model. For 4U 1626-67, 
M « 2.7 X 10^*^5/5 with hy^B^ ^ 2 y. lO^^g. We estimate M ~ 6 x lO^^g/s and 
hy^B^ ~ 5 X lO^^G for GX 1+4, and M ~ 1 x lO^^^/s and h^^B^ ~ 2 x lO^^^;. foj, 
OAO 1657-415. Reliable detection of QPOs before and after torque reversal could 
directly test the model. 



Subject headings: accretion, accretion disks — pulsars: general — stars: magnetic fields 
— X-rays: stars 



1. Introduction 

Sudden torque reversal events in some accretion-powered X-ray pulsars such as 4U 1626-67, GX 
1+4, and OAO 1657-415 have recently been detected (e.g. Chakrabarty et al. 1993,1996, 1997ab). 
The spin-up and spin-down rates are puzzlingly similar despite abrupt torque reversal. The 
torques remain largely steady before and after reversal, which plausibly indicates the existence of 
an ordered, stable accretion flow. These systems are distinguished from those showing random 
torque fluctuations seen in some wind-fed pulsar systems (e.g. Nagase 1989, Anzer k. Borner 
1995). In the torque-reversing systems, the mass accretion rate M appears gradually modulated 
with a small amplitude on a time scale ^ yr which is much longer than the typical reversal time 
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scale ^ days. The recently reported flux and spectral changes around the time of the reversal in 
4U 1626-67 shows that the M seems to change by about a few 10% in the X-ray emitting region 
close to the neutron star (Vaughn & Kitamoto 1997). 

Recently, Yi, Wheeler, & Vishniac (1997) |^ have suggested an explanation for the torque 
reversal phenomenon. The reversal is triggered by transition of the accretion flow from Keplerian 
rotation to sub-Keplerian rotation. The transition occurs when M crosses the critical rate 
Mc ~ 10^^ — 10^^5(/s. For magnetized pulsar systems, the inner region of accretion flow lies 
roughly at a radius i?o 5 x 10^5^/2^16^''''^"^' where Mig = M /IQ^^g/s and -B*,i2 is the stellar 
field strength in units of 10^^ G. This radius is similar to the size of a white dwarf, which strongly 
suggests that the transition could be similar to that of cataclysmic variables (e.g. Patterson & 
Raymond 1985, Narayan & Popham 1993). The details of the transition process remain unknown. 
The relevant time scale for the transition is likely to be the thermal time scale tth ~ (a^^i^)"^ or 
the viscous-thermal time scale t^t ~ {R/H)tth, where a is the conventional a viscosity coefficient, 
H is the thickness of the accretion disk at radius R from the star, and 0,k = {GM^, / R^Y^"^ is the 
Keplerian rotation rate for the stellar mass M*. The time scales tth < i-ut ^ lO^s for a ~ 0.3, 
R ~ 10^ cm, and Mig = 1 (Prank, King, & Raine 1992). After the transition, the rotation of the 
accretion flow becomes sub-Keplerian due to large internal pressure (cf. Narayan &: Yi 1995). 
The change of rotation induces an abrupt (^ day) decrease of torque exerted on the star by the 
accretion flow. The observed X-ray emission is not expected to be significantly affected since M 
change is small. 

We show how the proposed model tightly constrains the pulsar system parameters. The 
estimated parameters are largely consistent with other available estimates. We also derive a simple 
criterion which identifies candidate pulsar systems by their observable parameters. Combining 
this with an additional constraint from the quasi-periodic oscillation (QPO), we suggest a 
possible test of the model. We also discuss different types of torque variation and how they differ 
observationally from each other. The observed reversal events are likely to occur in pulsar systems 
near spin equilibrium with M ~ Mc ~ a few xW^^g/s. We take the neutron star moment of 
inertia = W^^g cm?, radius i?* = 10^ cm, and M^, = \AMq. The magnetic field is assumed to 
be a dipole configuration (Prank et al. 1992). The angular velocity of the star Vt^ = 2tx jP^. 

2. Accretion Flow Transition and Torque Reversal 

When the accretion flow rotation is Keplerian, the corotating stellar magnetic field lines 
interact with the rotating accretion flow with i^xiR) (e.g. Yi et al. 1997, Wang 1995 and 
references therein). In a steady state, the pitch of the azimuthally stretched magnetic field is given 



^An erratum is to be published in ApJL due to errors in the quoted parameters. The corrected parameters include 
substantially higher magnetic fields and mass accretion rates which are consistent with those in this work. 



-3- 



by 

B,{R) a nK{R) ' ^' ' 

where 7 < 1 measures the vertical velocity shear length scale. The dimensionless ratio bp = 7/0;, 
which determines the magnetic pitch, is not well constrained, but is likely to be order unity (e.g 
Wang 1995 and reference therein). The magnetically disrupted accretion flow gives an inner edge 
at R = Ro determined by 



S = 



Rc 



1 



RoV/'' 



R, 



-1 



(2-2) 



where Rc = {GM^P^ /4n'^y/^ is the Keplerian corotation radius and 

32 



~ (GM,)2/3 p7/3^ ^ 10 ^e//,lin,10 ^x,36- (2-3) 

We have defined B^jf n = bl^^ B^/IO^^G and = GM^M / R^. The torque exerted on the star is 

l-(8/7)(i?Ji?j3/2 
"^-6^° l-{Ro/Rcf'- ^'"'^ 

where Nq = M{GM^RoY/'^ (Yi et al. 1997, Wang 1995). The spin equilibrium = is achieved 
when Ro/Rc = Xeq = (7/8)^/'^. The pulsar spin evolution follows = —P^N/2TrI^ as M varies, 
which we model as a linear variation determined simply by dM/dt. If the accretion flow remains 
Keplerian with gradually varying M, the disruption radius and the resulting torque N varies 
according to eqs. (2-2), (2-4). If a pulsar system evolves from a spin-up (spin-down) state and M 
decreases (increases) gradually, it is possible that Ro/Rc = Xeq is reached at a certain M and then 
moves to Rq/Rc > (< Xcq) and hence spin-down (spin-up). As long as M variation is smooth 
and gradual the torque reversal event is expected to be smooth. Figure 1 shows an example of the 
smooth transition for 4U 1626-67, which is clearly distinguished from the observed abrupt reversal. 

When the decreasing M reaches Mc below which the accretion flow becomes sub-Kcplerian 
(Yi et al. 1997), the accretion flow-magnetic field interaction would change on a time scale ^ day. 
For sub-Keplerian rotation, fi/fi/^ = ^ < 1, the new corotation radius R'c = A^/^Rc and the new 
inner edge i?^ determined by 



5' 



7/2 r 
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(2-5) 



where 5' = 5 A ^he new torque on the star after the transition is 



7 l-(8/7)(i?-/i^-)3/2 
^ - 6^° 1 - {R'JR'cfl^ ' ^ ^ 

where N'q = AM^GM^R'^/^. As in the Keplerian fiow, N' = would occur when R'^/R'^ = Xeq- 
In a sub-Keplerian fiow supported by internal pressure, ^ = AUk is largely determined by the 
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ratio of magnetic to gas pressure in the accreted plasma. For equipartition between magnetic 

and gas pressures, with the pressure ratio, say, 0.05 — 1, we expect A ~ 0.14 — 0.4 for all 
a = 0.01 — 0.3 (e.g. Narayan & Yi 1995). For the Keplerian flow, the equilibrium spin period 
is Peq ~ 13{bp/10)^^'^ 31^12^16^"^ ^- -^^^ sub-Keplerian rotation, the new equilibrium spin 
period Pg^ = Peq/A > Peg. A system in a spin-up state with a restricted period ratio given by 
1 ^ P*/Peq ^ before transition, would evolve toward P^^ ^ P* after transition, i.e. a sudden 
torque reversal. Figure 1 shows a flt to the observed reversal event in 4U 1626-67 (cf. Yi et al. 
1997). Vaughn and Kitamoto (1997) report that M changes by ~ 20% before and after the the 
observed reversal. The model in Figure 1 is consistent with this finding. On the other hand, if a 
system has a pre-transition P^/P^g > A~^ that places it well into the spin-up regime, the transition 
could merely lead to a sudden decrease of spin-up rate, but not to spin-down (i.e. P* > Pg^), a 
behavior most likely in systems with smaller S*. Figure 1 shows a hypothetical pulsar similar to 
4U 1626-67 except with a smaller B^. Observations of this type of torque change could support 
the present accretion flow transition model. 

Quasi-periodic oscillation (QPO), if observed, could provide an additional constraint. Such a 
constraint would depend critically on the nature of the QPO mechanism. Adopting the widely 
used beat frequency model (e.g. Lamb et al. 1985), we get the QPO periods 



QPO 



= P* 



(p„/Pg)-3/2 - il \ p;,po = pJ{K/RT'^'-i\ ' (2-7) 



before and after transition, respectively. It remains unclear whether the conventional beat 
frequency model can explain the recently observed kHz QPOs in some X-ray pulsar systems. The 
QPO periods we use in the following section are much longer than than the kHz QPO periods. 
There exists no indication that the long period QPOs arise from other physical processes. 



3. Application to Observed X-ray Pulsar Systems 

The expressions derived so far can be readily applied to the observed systems. 

4U 1626-67 has P* 7.660s at reversal and the observed torques for the adopted /* 
are N 5.37 x 10^^gcm?/s'^ and N' —4.51 x 10^^gcm?/s^, or the observed torque ratio 
{N'/N)obs ^ -0.840 (Chakrabarty 1996, Chakrabarty et al. 1997a). Using the observed 

Pqpo ~ 25s during spin-up (Shinoda et al. 1990), eq. (2-7) gives Ro/Rc ~ 0.84 or 5 k, 2.3 from 
eq. (2-2). Using the observed N and the derived Rq/Rc in eq. (2-4), we get M ^2.7 X IQ^^g/s. 
Making use of the observed A^' and 5' = 5A~'^/^ , we solve eqs. (2-5), (2-6) to obtain A « 0.46, 
which suggests that the accreted plasma has an equipartition strength magnetic field. From eq. 
(2-3), we then find Pe// = tp^^P* ~ 1-7 x lO-'^^G. This field strength is smaller than the estimates 
P* ~ (6 - 8) X W^G obtained by Pravado et al. (1979) and Kii et al. (1986) based on the energy 
cutoff in X-ray spectrum and the energy dependence of the pulse shape. Although uncertainties 

1/2 

in the latter estimates of P* are not clear, a simple comparison suggests that bp is at most order 
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unity (cf. Wang 1995). The observed 0.7-60 keV flux ^ 2.4 x 10~^erg/s/cm? (Pravado et al. 
1979) gives the distance estimate (Lx/4:7rFx)^^'^ ~ 4.2 kpc, which is consistent with the previous 
estimates (e.g. Chakrabarty 1996). Eq. (2-7) suggests a possible P'qpq ^ 120s after transition. 
Detection of such a QPO could directly test the model. 

GX 1+4 reversed torque from spin-down to spin-up around the spin period ~ 122.15s 
(Chakrabarty 1995, Chakrabarty et al. 1997b) with the measured torques N ~ 3.77 x 10^^ (/cm^/s^, 
and N' Ki -3.14 x 10^^ gcm'^ / s'^ , which give iN'/N)obs ~ -0.833. A QPO of period Pqpo ~ 250s 
was reported during the 1976 spin-up (Doty et al. 1981), but none near the recent reversal 
episode. The applicability of this QPO is therefore questionable. Nevertheless, assuming that the 
1976 spin-up state is similar to the recent spin-up, we can adopt this QPO to make estimates 
for GX 1+4. The spin period around 1976 is close to the recent spin period, which may help 
to justify our use of the 1976 QPO data. Following the same procedure as for 4U 1626-67, we 
estimate that Ro/Rc ~ 0.77, S ~ 1.2, and M ~ 2.9 x lO^^g/s. We get A ~ 0.12, which again 
indicates a magnetic field strength near equipartition. Finally from eq. (2-3), -Bg// ~ 3 x lO^^G or 
5* ~ lO^^G for bp ~ 10. The observed flux 2 x lO"^^ erg/s/cm^ in the range of 20-60 keV 

(Chakrabarty et al. 1997b) gives a distance estimate of ~ Ibkpc for the derived M whereas for 
the Doty et al. (1981) value oi Fx ^ 8 x W^^erg/ s/cm^ in the 1.5-55 kcV range gives a distance 
estimate of ~ 2.4fcpc. During the recent spin-down, Pqpq ~ 6160s would be predicted. This 
prediction, however, is much less certain than that for 4U 1626-67 due to the lack of reliable Pqpo 
measurement near the recent spin-up. We note that the recent GX 1+4 reversal is considerably 
more gradual than the 4U 1626-67 event. It is possible that the smooth transition type shown in 
Figure 1 could be relevant in this case (Yi et al. 1997). 

OAO 1657-415 has an observed pulse period 37.665s at the time of reversal. Recent 

observations give 

AT ~ 4.40 X l(y^^gcm?/s^ and A^' ^ -1.06 x lOr"^ gcni^ / (Chakrabarty et 
al. 1993) or {N' /N)obs ~ —0.241. There exists no reported QPO for this system. If we take 
A = 0.4 based on our estimates in 4U 1626-67 and GX 1+4, the observed torque ratio {N' /N)oi,s 
corresponds to (5 ~ 1.2 or Rq/Rc ~ 0.77 as in GX 1+4, which gives M ~ 1.0 x 10^'^g/s. Therefore, 
we get Beff ~ 2 X lO^^G or ~ lO^^G for bp ~ 10. For P* ~ 37s, Ro/Rc ~ 0.77, and A = 0.4, 
we expect Pqpo ~ '75s and Pqpo ~ 430s before and after reversal, respectively. For a detected 
flux of ~ W~^erg/s/cm'^ (White & Pravado 1979, Mereghetti et al. 1991) and the estimated 
M ~ 1 X 10^'^ g/s, the distance ~ 12kpc. 

We caution that the predicted QPO frequencies may not be easily detectable. During 
spin-down, the accretion flow thickens considerably and the pulsed polar emission may be 
significantly diluted by scattering (Yi et al. 1997). The long period QPOs may appear as 
occasional flares with repetition time scales of ^a few 10^ s. 
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4. Necessary Condition for Torque Reversal in Systems Near Spin Equilibrium 



The parameters derived above can be understood by a simple expression that constrains the 
pulsar parameters and essentially determines the pulsar magnetic field. We assume that the system 
reverses its torque from spin-up to spin-down. Just before reversal, we take Rq/Rc = x^q — e with 
e < Xeq- Expanding eqs. (2-2), (2-4) to first order in e/xgq and combining we get 



N 



7No e/x. 



eq 



4 1-4' n7-4x7q 



3/2 X 



1 ~ x] 



3/2 
eq 



7/2 



(4-1) 



After transition to the sub-Keplerian state, we take R'^/R[ 
expand eqs. (2-5), (2-6) to first order in e' /xeq and combine them to derive 



Xeq + e' also with e' < x^q- We 



N' 



7N' 



3/2 



7/2 
Xeq 



5'-l 



(4-2) 



4 l-x%''^eq 2(7-4x^f) 
The ratio of the sub-Keplerian spin-down torque {N') to the Keplerian spin-up torque (N) is then 

g{Xeq)6' - 1 



N 



I -I- 9{Xeq){^ Xeq ) _ 



4(7/4 - x%') 



1 - g{Xeq)S 



(4-3) 



where g{xeq) = (1 — x'jq^jxjq. We have made use of the first order expression 



^=^4/3 [i + + ^^4/3 



'^eq 



X. 



eq, 



4(7/4 - xli') 



For Xeq = (7/8)2/3, we get 



N 



1 + 6.10 X 10-^(6' - 6) 



0.171(5' 



0.171(5 



(4-4) 



(4-5) 



Given the observational constraint N' ^ N, eq. (4-5) indicates that small values of ^ ~ 0.1 
are unlikely, which is consistent with the result A > 0.3 in the previous section. This in turn 
suggests that the magnetic to gas pressure ratio is roughly at the level of equipartition (cf. 
Narayan &; Yi 1995). Prom eqs. (4-1), (4-2), the first necessary condition for the torque reversal is 
N'/N < or 



^7/3 



,7/2 



<6< 



7/2 

'^eq 



, 3/2 " " " , 3/2 
-L Xeq J- Xeq 



or for Xeq = (7/8)2/3 



^„,„ = 5.86A^/3 <6 <5I 



(4-6) 



(4-7) 



The second condition comes from the observational requirement, |iV'/A^|^j^ < 1, which gives (cf. 
eq. 4-5) 

1 + ^4/3 

(5 < 5.86- (4-8) 

1 + A-^ ^ ' 
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Finally, the observed continuous X-ray emission throughout the whole reversal process requires 
the accretion to be continuous. This translates into the condition Rq/Rc < 1 and R'^/R'^ < 1 in 
order to avoid an angular momentum barrier at the corotation radius. By expressing e' in terms 
of 5, we get 

5 < 13.39^^/=^ = Smax- (4-9) 

For A < 0.51, 5 

max replaces the upper bounds in eqs. (4-7), (4-8). For A > 0.51, the upper bound 
on S is given by eq. (4-8) as shown in Figure 2. As A — 1, the allowed parameter space shrinks 
rapidly, which simply indicates that the rotation gets closer to the Keplerian and the spin-down 
becomes more difficult when the mass accretion rate is fixed. Therefore, the rotation needs to be 

substantially sub-Kcplcrian after transition in order for the reversal to occur. For most of the A 
values, 6's range is limited to be within a factor of ^ 2. Therefore the condition for torque reversal 
becomes 0,711^ ^max and the magnetic field should be close to 

B^ff ^ 7 X W'Hy'Ll%P:;,lG (4-10) 

where dlll^ = 2.4^^/6 < 5V2 < 3.7^7/6 = 

Srnax- If the accrcted plasma is magnetized with 
near equipartition between magnetic pressure and gas pressure, the sub-Keplerian flow with 
^ ~ 0.1 — 0.4 (Narayan &; Yi 1995) suggests that (5^/^ is essentially a constant of order unity (cf. 
Figure 2). Therefore, eq. (4-10) provides an estimate of the magnetic field within a factor of ~ 2. 

If the distance to the pulsar and an independent estimate on (say through X-ray spectrum and 
emission) are known, this relation can constrain bp = (i?e///i?*)^- The pulsar parameters derived 
in the previous section are consistent with eq. (4-10). This can be checked as follows. Using eqs. 
(4-5), we derive 

0.1707 A-i - {N>/NU, ■ ^' ' 

For a given A, we confirm that the observed torque ratio gives 5's which are close enough to the 
derived values in the previous section. Therefore we conclude that eq. (4-10) can be reliably used 
to identify the characteristics of the accreting pulsar systems that should show torque reversals. 

The condition for the torque reversal, 5.86^4^/'^ < 6 < 13.39A^/^, for A < 0.51 allows 
a very tightly constrained physical parameters of the pulsar systems. Based on the results 
in the previous section, our model suggests A « 0.4 and therefore 0.69 < 5 < 1.58. Then, 
in eq. (4-10), 0.8 < 5^/^ < 1.3 is determined accurately. If this is the case, we can derive 
Beff ~ 7 X 10^^ Llf^QPy^QG essentially without uncertainties in S. Since S oc A'^^^, B^ff oc A'^^^ 
has a crucial dependence on A. 



5. Discussion 

The current analysis and the study of Yi et al. (1997) suggest that the observed torque 
reversals occur when the accretion flow makes a transition at a critical rate Mc ~ 10^^ — lO^^^/s. 
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This points to an interesting connection between X-ray pulsar systems and other compact 

accretion systems such as cataclysmic variables and black hole transients. We speculate that the 
spectral transitions seen in the latter systems may be due to a common accretion flow transition 
mechanism (cf. Narayan & Yi 1995). The observed torque reversal could be a signature of a pulsar 
system near spin equilibrium with M gradually varying near Mc- The proposed mechanism tightly 
constrains the pulsar field strength. If a pulsar does not have an appropriate field strength (eq. 
4-10), it would not undergo a torque reversal even if M passes through Mc causing the transition. 
If all systems have a similar (Narayan & Yi 1995), the pulsar magnetic field strength can be 
determined despite a major uncertainty in the distance. It is possible that some systems with 
P* > Peg abruptly change their torques without reversing sign at M ~ Mc- These types of sudden 
torque change are distinguishable from the smooth torque variation in a Keplerian accretion flow. 
If QPOs can be detected near torque reversal, the proposed model could be further tested in 
detail. 

Although the spectral and flux changes are generally small in the torque reversing systems 
(Chakrabarty 1996), which indicates the small changes in M, at least in 4U 1626-67, some spectral 
and flux changes have been reported (Vaughn & Kitamoto 1997). The hardening of the spectrum 
after the transition to spin-down could be attributed to the accretion flow transition itself (Yi 
et al. 1997). It is however unclear if the hot accretion flow with a large vertical thickness could 
directly affect the X-ray spectrum through scattering since the expected scattering depth through 
the hot flow is at most a few xlO~^ for M ~ 2 x lO^^g/s (Narayan & Yi 1995) appropriate for 4U 
1626-67. 

Gradual M modulation on a time scale ranging from ~ yr to a few decades is required for 
the proposed torque reversal. Such a modulation can be driven by several known mechanisms, (i) 
In the case of GX 1-1-4, the binary orbital motion on a time scale ~ yr is plausible (Chakrabarty 
1996), which could cause M variation through orbital modulation of mass transfer. In 4U 1626-67, 
an orbital time scale longer than a few years is unlikely (Rappaport et al. 1977, Joss et al. 
1978, Shinoda et al. 1990). OAO 1657-415 has an orbital time scale of ~ lOd (Chakrabarty 
et al. 1993). (ii) Several precession time scales involving the pulsar and the secondary remain 
viable (e.g. Thorne et al. 1986). (iii) X-ray irradiation- induced mass flow oscillation (Meyer 
& Meyer-Hofmeister 1990) could provide long time scale (> month) oscillations when certain 
conditions such as disk size and a are met. Interestingly, in 4U 1626-67 the observed optical 
pulsation frequency is the same as that in X-rays, which has been attributed to the reprocessing of 
X-rays by the accretion disk (Ilovaisky et al. 1978, Chester 1979). (iv) The disk instability model 
(e.g. Smak 1984) has not been applied to neutron star systems and long term M variation due to 
the disk instability needs further investigation. 

We acknowledge receiving useful references from L. Bildstcn and D. Chakrabarty. JCW 
appreciates the hospitality of the Institute for Theoretical Physics which is supported by NSF 
grant No. PHY94-07194. This work has been supported in part by SUAM Foundation (lY) and 



- 9 - 



by NSF grant AST95-28110 (JCW). 



REFERENCES 

Anzer, U. & Borner, G. 1995, A&A, 299, 62 

Chakrabarty, D. 1996, Ph.D. Thesis, California Institute of Technology 
Chakrabarty, D. et al. 1993, ApJ, 403, L33 
Chakrabarty, D. et al. 1997a, ApJ, 474, 414 
Chakrabarty, D. et al. 1997b, ApJ, 481, LlOl 
Chester, T. J. 1979, ApJ, 227, 569 

Doty, J. P., Hoffman, J. A., & Lewin, W. H. G. 1981, ApJ, 243, 257 

Prank, J., King, A. R., & Raine, D. 1992, Accretion Power in Astrophysics (Cambridge: Cambridge 
University Press) 

Ilovaisky, S. A., Motch, C, & Chevalier, C. 1978, A&A, 70, L19 

Joss, P. C, Avni, Y., & Rappaport, S. 1978, ApJ, 221, 645 

Kii, T., Hayakawa, S., Nagase, F., Ikegami, T., & Kawai, N. 1986, PASJ, 38, 751 

Lamb, F. K., Shibazaki, N., Alpar, M. A., & Shaham, J. 1985, Nature, 316, 239 

Mereghetti, S. et al. 1991, ApJ, 366, L23 

Meyer, F. & Meyer-Hofmeister, E. 1990, A&A, 239, 214 

Nagase, F. 1989, PASJ, 41, 1 

Narayan, R. & Popham, R. 1993, Nature, 362, 820 
Narayan, R. & Yi, I. 1995, ApJ, 452, 710 
Patterson, J. k Raymond, J. C. 1985, ApJ, 292, 535 
Pravado, S. H. et al. 1979, ApJ, 231, 912 
Rappaport, S. et al. 1977, ApJ, 217, L29 
Shinoda, K. et al. 1990, PASJ, 42, L27 
Smak, J. 1984, PASP, 96, 5 

Thorne, K. S., Price, R. H., & MacDonald, D. A. 1986, Black Holes: The Membrane Paradigm 
(New Haven: Yale University Press) 



Vaughn, B. A. & Kitamoto, S. 1997, preprint (|astro-ph/9707105D 

Wang, Y.-M. 1995, ApJ, 449, L153 

White, N. E. & Pravado, S. H. 1979, ApJ, 233, L121 



- 10- 

Yi, I., Wheeler, J. C, & Vishniac, E. T. 1997, ApJ, 481, L51 (erratum to appear in ApJL) 



This preprint was prepared with the AAS lAT^^X macros v4.0. 



Figure Caption 



Figure 1: Left columns: Torque reversal event of 4U 1626-67. The solid line is the 
combined BATSE and pre-BATSE spin period history adopted from Chakrabarty (1996). The 
dotted lines correspond to the model with the accretion flow transition where = 5.4 x 10^^ G, 
the mass accretion rate decreases from an initial rate of Mj = 3.2 x 10^^ g/s at a rate 
dM/dt = — 6 X IQ^'^g/s/yr. The accretion flow makes a transition at Mc = 2.4 x IQ^^g/s. 
The dashed line corresponds to the Keplerian model without accretion flow transition where 
= lO^^G, Mi = 5.9 X W^^g/s, and dM/dt = -2.35 x lO^^g/s/yr. The observed event is 
noticeably more sudden than the latter model. Right Column: A hypothetical pulsar similar to 4U 
1626-67 (with abrupt reversal) except i?* = 3 x 10^^ G. When > Peq/A (i.e. well into spin-up), 
the accretion flow transition results in abrupt decrease of torque without reversal (dotted line). 
The dashed line is the corresponding Keplerian model without transition. The parameters are 
bp = 10 and A = 0.46 (see text). 



Figure 2: Allowed region in the S — A parameter space for the torque reversal. The region 
between the two lines are allowed for the torque reversal. 
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